from IPython import display
display.Image('team.png', width=700)
display.Image('./old_work_flow.jpg')
# data: http://www.amber.utah.edu/AMBER-workshop/London-2015/DNA-tutorial/
import pytraj as pt
traj0 = pt.load('md.nc', 'dna.prmtop')
traj0
traj = traj0.autoimage()['!:WAT']
traj
# compute rmsd and convert raw data to pandas' DataFrame
data = pt.rmsd(traj, ref=0, mask=':1-14&!@H=', dtype='dataframe')
data.head(5)
%matplotlib inline
data.plot()
data.hist()
traj = pt.iterload('md.nc', 'dna.prmtop')
data = pt.pmap(pt.rmsd, traj, ref=0, mask=':1-14&!@H=', n_cores=8)
# serial: data = pt.rmsd(traj, ref=0, mask=':1-14&!@H=')
from IPython import display
display.Image('bench_pmap_casegroup.png', width=500)
traj2 = pt.iterload('tz2.nc', 'tz2.parm7')
energies = pt.energy_decomposition(traj2, igb=8, dtype='dataframe')
energies[['bond', 'angle', 'dihedral', 'gb']].head()
# XTC
import mdtraj as md
t0 = md.load('monolayer.xtc', top='monolayer.pdb')
coordinates = t0.xyz.astype('f8')
traj = pt.Trajectory(xyz=coordinates, top='monolayer.pdb')
pt.center_of_mass(traj)
help(pt.calc_center_of_mass)
display.Image('./3pqr_nglview.png')
display.Image('./trpzip2_nglview.png')